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Abstract. The 'ultra-discrete limit' has provided a link between integrable 
difference equations and cellular automata displaying soliton like solutions. In 
particular, this procedure generally turns strictly positive solutions of algebraic 
difference equations with positive coefficients into corresponding solutions to 
equations involving the "Max" operator. Although it certainly is the case that 
dropping these positivity conditions creates potential difficulties, it is still pos- 
sible for solutions to persist under the ultra-discrete limit even in their ab- 
sence. To recognize when this will occur, one must consider whether a certain 
expression, involving a measure of the rates of convergence of different terms 
in the difference equation and their coefficients, is equal to zero. Applications 
discussed include the solution of elementary ordinary difference equations, a 
discretization of the Hirota Bilinear Difference Equation and the identification 
of integrals of motion for ultra-discrete equations. 



1. Introduction 

Integrable nonlinear partial differential equations are of special interest in part 
because of their 'soliton solutions' which exhibit particle-like behavior [1,3,7,23]. 
The link between these equations and the cellular automata that exhibit superfi- 
cially similar solutions [2,11,13,18] was the observation in [20] that solutions to the 
discretizations of the soliton equations become solutions to these cellular automata 
under the 'ultra-discrete limit'. Loosely speaking, this works because of the fact 
that for a, 6 G R 

(1) lim e In (e""/^ + e''/^) = Max(a, b). 

This formula can be understood intuitively by noting that the term corresponding to 
the smaller of the two exponents will be relatively insignificant when e is very small 
and that with it eliminated the entire expression is simply equal to the numerator 
of the other exponent. It can also be easily verified using L'Hopital's rule. 

In applying to difference equations, previous authors have been very cautious 
allowing negative coefficients in front of the exponentials, and have not been con- 
cerned at all about the question of how allowing a and b in the formula above to 
depend on e might affect the outcome. In fact, they are correct that letting those 
coefficients be negative may invalidate the prediction that the resulting limit will 
merely be given by the "Max" . In this paper, we will be especially interested in the 
case in which the linear combination of exponential functions does involve negative 
coefficients. Our viewpoint is that it is worthwhile to investigate this situation to 
understand what can go wrong. 

Theorem 12.11 demonstrates two surprising facts: (a) Even if negative coefficients 
are involved, a formula analogous to Q continues to work "most of the time" . 



1 



Negativity Not Problem For Ultra-discrete Limit 



A. Kasman/S. Lafortune 



(b) In determining when it might fail it is necessary to explicitly consider the rate 
of convergence of the exponents under the limit in e. This theorem is illustrated 
by an application to two simple, ordinary difference equations in Section |3| In a 
more sophisticated application, the theorem is used to guarantee the production of 
a solution to an ultra-discrete limit of the Hirota Bilinear Difference Equation in 
Section 0] Finally, these results are applied to the question of when integrals of 
motion for difference equations survive the ultra-discrete limit in Sectional 

Both to refer readers to other work that may be of interest and to emphasize 
the differences with this paper, we would like to conclude this introduction with 
a discussion of three other papers. Equation |^ has been previously generalized 
in [21] to the case in which a and h are complex constants. This is of great interest, 
but clearly quite different from what we do below by allowing the argument of the 
logarithm to be an arbitrary linear combination of exponentials and letting the 
exponents be real functions of e. Also, [8] presents work related to the goals of 
the present paper by addressing the question of how one can take an ultra-discrete 
limit of a solution to a difference equation which takes on negative values. Their 
approach involves an alternative method of introducing dependence on e into the 
solution so that the new variable can be real even in the case that the original 
variable was negative. However, in doing so they are careful to be certain that 
they always apply ^ only in the case that the coefficients of the exponentials are 
equal to one and make no explicit mention of the possibility that the results may 
be affected by the dependence of the exponents on e. So, that paper takes what 
might be viewed as an "orthogonal" approach to the question to the one that will be 
pursued below. It might be of interest to attempt to combine the two, by applying 
Theorem 12. II to the s-ultra-discrete limit as defined in that paper. Finally, perhaps 
most closely related to this paper in spirit, [12] generalizes the procedure of ultra- 
discretization by formally extending the Max-Plus algebra through the introduction 
of an extra Z2 component. This extension of the Max-Plus semi-ring to a full ring, 
as is normally done in the Max-Plus community by moving to the symmetrized 
ring Smax [4], means that the resulting equations are not ordinary ultra-discrete 
equations (the Max operator and even equality have slightly different meanings). 
In contrast, in this paper we address the question of when negativity fails to create 
any problems under the original definition of the ultra-discrete limit and integrable 
ultra-discrete equations. 

1.1. Notation. In the following, a superscript of e will indicate that the parameter 
depends dijferentiably on the real variable e in some open interval interval (0,/3). 
Thus would be an e dependent bi-infinite sequence and could be a solution to a 
difference equation in which the coefficients also depend on e, and we would only 
call it a solution if the equation were true not only for all n but also for all e € (0, /?). 
(E.g. — en satisfies = (ufj^ — e^). 

1.2. An Enlightening Example. As we will see in Theorem 12.11 below, it is 
possible to generalize ^ to the case in which the argument of the logarithm is a 
more general linear combination of exponentials whose exponents depend on e. As 
one might guess, the value of the limit turns out to be equal to the maximum of 
the limits of the numerators of the exponents. However, the theorem can only be 
applied subject to the condition that a certain expression involving the coefficients 
and the rates of convergence is not equal to zero. Following is a cautionary example 
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whose purpose is to illustrate the fact that failing to consider this condition is liable 
to produce incorrect results. 
If we consider the expression 

r = e"i/^ + 7e"2/= 

then using the previously stated argument that the term with the lower power can 
simply be ignored, we may naively think that 

(2) lim elnr = Maxfui, U2) where ut = lim u^. 

However, ^ may be false for some choices of uf and 7. 
For instance, if we let 

u\ = a + 2£\n {e^''/^ + + e-2fe/£^ i4 = a + eln2 and 7^-4 
where fc > and a is arbitrary, then 

ui = U2 = a. 

However, despite the predictions of equation Q, one can check that in fact 

lim elnr = a — k. 

The total avoidance of negativity in the literature on "ultra-discrete limits" might 
suggest that it is the fact that 7 < that is causing a problem here, but that would 
be a drastic oversimplification. The fact is that for any choice of 7 ^ —4 the 
prediction of jSjl would have been correct and the limit would have the value a as 
expected. So, an explanation of this phenomenon would have to shed light on the 
question of why this occurs only at 7 = —4. (Such an explanation is provided in 
the following section.) 

This example is not as artificial as it may at first seem since it will arise in 
a more general form as part of an investigation of a simple ordinary difference 
equation in Section 13 It demonstrates that the coefficients in the linear combina- 
tion of exponentials can have unexpected effects on the "ultra-discrete limit". It 
also demonstrates that the natural generalization of provided by (0) is correct 
generically in the sense that it is only for rather specific choices of coefficient that 
it fails. 



2. A More General Limit 

This section states and proves a generalization of which will later be applied 
to questions involving difference equations. We wish to consider instead the case in 
which the argument of the logarithm, r, may be an arbitrary linear combination of 
exponentials and in which the numerator of the exponents can also depend upon 
e. The following theorem gives conditions that are sufficient (but not necessary, as 
we will see in Section I5!T|) to conclude that the limit of elnr is still equal to the 
largest limit of the numerators of the exponents. 

Theorem 2.1. Let 

n ^ 

(3) r = ^7,e^ 
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where 7, e C are arbitrary, non-zero constants and wf are real valued functions 
such that the limits 

Mi = lim i = 1, . . . ,n 

£^0+ ' 

are finite. Let 

m = Max(Mi, . . . , M„) and 3 = {i\Mi = OJl} 

denote the largest value of these limits and the set of indices for which this maximum 
is achieved, respectively. 
Then 

lim e In r = 9Jl 

if either 

i. \3\ = 1 (the maximum is achieved only once) 
or 

ii. the limits 

Hi = lim^(w;f - 9Jl)/e for ie3 



are finite and D ^0 where 



D = > 7.e'^^ 



ie3 



lim elnr = lim eln 

£^0+ £^0+ 



Proof In the case \3\ = 1 let J = {io} so that Mj^ = 9Jl > for i 7^ io. Then 
= £lne^/^+eln^7ie(-.-^)/^ 

(then since wf — SOt < for i ^ zq and e sufficiently small) 

= lim £ln7jo = SDT. 

£— >0+ 

The situation is more complicated in the case that \3\ > 1. We now assume 
that D 7^ and show that again the limit is simply equal to the maximum of the 

exponents. 

First, we wish to show that D ^ implies that r is non-zero on a small open 
interval (0,/3) and so it makes sense to consider the function elnr. (There is no 
need to require t > 0, however, since the proof to follow continues to be valid in the 
case that this function takes on complex values or even if it becomes multivalued 
as a consequence of analytic continuation.) 

D = ^7,e''* = lim^^7,e«-^)/^ 

(then since — 971 < for i J and e sufficiently small) 

n 

= lim yjj<-^y^= lim e-™/^T 

£^0+ £^0+ 

Suppose £i > is a sequence that converges to and with the property that t = 
for e = ej. Then, since the limit above converges to D, it must be the case that 
D = 0. Since D ^0, no such sequence of {si} can exist. 
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Now, we simply compute that 



lim e In r ~ lim e In 

6^0+ £-»0 + 



n 

lim elne™/^+ lim eln V e^""'"™)/" 



i=l 



= 971 + lim e\n\Y. e^"'' + E e^"'''"^^/" 
= 971+ lim e In (£> + 0) = 971. 



□ 



Note that, apart from the restriction that _D 7^ 0, the value of the limit is 
independent of the choice of coefficients 7;. We will take advantage of this feature in 
Sections 13^ El to effectively apply the procedure of ultra-discrete limits to integrable 
difference equations even in the absence of any positivity restrictions. Furthermore, 
note that as a consequence of this lack of dependence on the value of the coefficients, 
the limit will be unaffected by an extra common factor of —1. In particular, we 
have the following corollary: 

Corollary 2.2. If the conditions of Theorem \2. 1\ are met then 

lim elnlrl = lim elnr. 

Again, this ability for us to take the absolute value prior to taking the ultra- 
discrete limit will be useful in our attempt to apply this result in the case of differ- 
ence equations. 

2.1. Continuation of Earlier Example. It is important to note that if neither 
conditions (i) or (ii) are met, then the value of the limit cannot necessarily be 
determined through knowledge of the limits Mi of the functions in the exponents. 
This is demonstrated by the example from Section 11.21 in which the value of the 
limit of the entire expression, a — k, depends on the parameter k not appearing in 
the limits of the functions in the exponents when 7 — —4. 

We can now apply Theorem 12.11 to explain the significance of 7 = —4 in that 
example. Again let 

T = e"i/''-|-7e"^/'^ u\ = a+2e In [e.-^'" + VsTT^^fcT?) and = a+e In 2 

but consider 7 arbitrary. The paramters /i^ measuring convergence are then 
/ii — lim iu\ — a) / £ — \nS and ^2 = lim (uS ~ a) le = In 2. 

Hence, we find that 

i? = e'"«H-7e'"2 = 2(4-f 7). 
Since this is only equal to zero at 7 = —4, the "mystery" is resolved. 

2.2. An Example vi^here D Does Not Exist. To illustrate the role of condition 
(i) in Theorem l2.1l we briefly point out that for 

^ ^ e"i/'^ - e"2/^ uf=esin(l/e) u| = -2 + 3e 

one has 971 = lim£^Q+ u\ — Q > limg^Q+ U2 = —2. Since |3| = 1 we know that 
lim£_,o+ ehiT = even though the limit defining /ii has no value. 
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3. Application: Ordinary Difference Equations 
Consider the difference equation 



N I d \ N I d 



(4) E n«+.)'"M =E «^ n^^' 

We will generally be interested in the case niij , rhij G Z for the sake of preserving 
the integrality of solutions under the ultra-discrete limit. 

Definition 3.1. We say that a solution to ^ is a candidate for the ultra- 
discrete limit if: 

• ^n-d^ • ■ • ' ^n+d' '^l ' '^'^^ ^'^^^ non-zcro real values for e G (0, /3) for some 
small, positive number (3; 

• The limits Un = lim eln|a;5, 1, Ai = lim elnlafl and Ai — lim elnldil 

£^0+ £^0+ e^0+ 

exist and are finite; 

• The limits 

{e\n\al\+eJ2j=-dTT^^j^^\<+j\) , . ^ 
fii = lim for I £ J 

£-»0+ £ 

and 

(eln|af| +e^^^_^myln|x^_^-|) -OT . . 
jli = lim for t E J 

exist and are finite where 

N ^ 

m = Max{A, + 2^ m,jUn+j) 

J = -d 

and 

d 

1 = {i\A^ + ^ rriijUn+j = Tl} 

j=-d 

is the set of indices for which this maximum is attained, and similarly, 

S!Jl = Max.{Ai + rhijUn+j) and 3 — {i\Ai + rhijUn+j = Tl}. 
j=~d j=-d 

Theorem 3.2. Suppose is a solution to the algebraic difference equation (0J) 
which is a candidate for the ultra- discrete limit (see Definition \3.1\ for terminology 
and notation). The ultra- discrete limit of xf-^^ 



lim £ In la;!: 



is a solution to the equation 



N jL N 



(5) Max(A, -f ) = Max(i, + Y m^jUn+j) 

j=-d j=-d 

if the maximum value on each side is never attained simultaneously for more than 
one value of the index, or as long as both 

(6) J2 '^^^^^ ^0 ""'^ E ^ ° 
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where 



are the signs of the corresponding term in the equation. 
Proof. Letting 



„ £in c4 n «+,)"'^- 

j=-d 

we can write the leftside of ^ in the form 



N 



If PI = 1 or Z) ^ 0, then Theorem 12 . II tells us that the limit of elnr will be 

N 

Max( lim wf). 

But that D is non-zero is exactly the first of the two requirements in ((^J) and so 



lim wf 

e-*0+ ' 



= lim I £ In a7 

£^0 + 



e ^ myln<_^j- 



— Ai ~\- ^ ^ rriijUn+j 

J = -d 

and that |3| — 1 is equivalent requiring that the maximum on the Icfthand side of 
(|SJl is achieved for a unique value of the index. 

Similarly, we conclude that the limit of e times the log of the righthand side of 
the equation is Max(4i -I- ''^I'ij'^n+j)- If 2;^ is a solution then these two sides 

are always equal and hence their limits are equal as well. □ 

3.1. Basic Example. Consider the difference equation 



(7) 



Let = 3 because there are three terms on the left, N — \ because there is one 
monomial on the right, and d = 1 because we only see n shifted up and down by 
this much. The monomials on the left are given by letting 



CTl = 1, r71l,_l ~ TOl o ~ 0. TOl l = 1, 



and 



0-3 = -1, 7713 1 = 1, W3 n = TO3 1 = 



The monomial on the right is given by letting 

(7i = 1, ol\ = 2e^^^ , TTii.-i — mi l = and mi.o = 1. 
A solution to this equation is given by 
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since expanding the lefthand side of jT)) with this definition for gives an expres- 
sion equal to the righthand side for any n or e. 

Trivially, this same solution also solves the difference equation 

(8) e^''^^x'^_i — e'^'^^'^x'^_i — 2e''^^x^ — x^^i- 

In fact, in most circumstances one would want to consider this equation to be 
completely equivalent to 10 since one can simply add a;^_|_j to both sides of one 
to make it into the other. However, for the purposes of this paper they are quite 
different. 

Regardless of whether we are working with the difference equation in the form 
(0) or lO, the solution xf^ = (e^/"^ + e'^/'^) is the same. Then our solution m„ to 
the ultra-discrete equation should take the form 

M„= lim £ln(e''/^ + e'=/^)" = n lim £ ln(e^/^ + e^^/^) = n Max(6, c). 

Ignoring condition (|SJ| of Theorem 13. 21 we get from that u„ should satisfy 

(9) Max(M„+i, 26 4- u„-i, 2c -I- Un-i) =b + Un 
and from JHJ we get that it should satisfy 

(10) Max(26 + u„_i, 2c -I- Un~i) = Max(6 -f- m„, u„+i). 

However, it is not the case that m„ satisfies each of these equations regardless of 
the choice of b and c. In particular, condition 10 may not be met and, as we will 
see, this causes u„ not to solve (0 when b < c. 

Look at the arguments of Max on the left hand side of and consider the 
question of which of them are equal to the maximum value. If 6 > c then the terms 
have the values (n -I- 1)5, (n -t- 1)6 and (n — 1)6 -I- 2c respectively so that the first 
two terms are equal to the maximum value. (Using the notation of the theorem, 
Tl = (n+ 1)6 and 3=1,2.) All of the /ii are zero (which means that these solutions 
converge quickly to their limits) and so all of the exponential terms are just equal 
to one. Now, we must look at the sum of ci and (72 to make certain that it is not 
equal to zero. Since cti + (T2 = 2 7^ the condition is met. And, in fact, one can 
easily verify that if 6 > c then m„ = nb is indeed a solution to 0. However, when 
6 < c then things are quite different. In that case, VJl = {n + l)c and 3 = {1, 3}, but 
(Ti + (73 = and so Theorem 13.21 does not predict that u„ will solve the equation. 
In fact, it does not because then the lefthand side has the value {n+ l)c while the 
rightside has the value nc + 6. 

In contrast, the theorem predicts that u„ will solve H1U|I when b > c and when 
b < c. In particular, if 6 > c then 971 = Max(26 + Un-i, 2c + u,i-i) — {n + 1)6 and 
3 = {1} while if c > 6 then dJl = {n + l)c and J = {2}. Thus, we conclude from 
the theorem that u„ should satisfy H1Q(I as long as 6 7^ c. In fact, it can be checked 
explicitly that u„ is a solution to (|10|l in the case b=c as well. This is true even 
though the conditions of the theorem fail to be satisfied in that case. This shows 
that the conditions of the theorem, while sufficient, are not necessary. 



3.2. Another example: Demonstrating the Role of Rate of Convergence. 

Consider the difference equation 
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Considering the survival of solutions of this equation under the ultra-discrete limit 
illustrates the importance of the parameters which measure the rate of conver- 
gence of the terms. 

This has a solution of the form 



Consequently, we would expect 



lim elnx'L 



nA A > 
A < 



to solve 

Max(u„+i, u„_i) ^ X + u„. 

And indeed, if A > this is true. However, it fails when A < 0. This could be 
predicted by the Theorem. Condition JHl is trivially met when A > because the 
exponents on the left do not have the same limit. However, when A < then both 
exponents go to and /ii = ln2 and hence the sum is = e'"^ — e'"^ = 0. 

4. Application: Hirota Bilinear Difference Equation 

The Hirota Bilinear Difference Equation (HBDE) [10,17,22] for a function x^{a, 6, c, 
of four variables is 



(A4 - A|)(A2 
+ (A|-AD(A| 
= (A|-A|)(A| 



Xl)x%a, b,c+l,d+ l)x''{a + 1, 6 + 1, c, d) 
A|)x^(a + l,b,c,d+ l)x%a, b+l,c+l,d) 
Xl)x%a + l,b,c+l, d)x''{a, b+l,c,d+ 1). 



As a consequence of the main result in [5], a solution to this equation can be 
constructed as 

x'ia, b, c, d) = det (P • (/ - Af S*)" ■ (/ - A^S')^ • (/ - X^Sy ■ {I ~ XlSf ■ n) 

where is any n x k matrix (fc < n), 5 is any n x n matrix such that the lower-left 
k X (n — k) block has rank one, / is the n x n identity matrix, and P = (0 Ik) is 
the k X n matrix which has the k x k identity matrix as a block and all zeroes to 
the left of it. 

It is interesting to note that even though x^ {a,b,c,d) has dependence upon e 
inheritted from the presence of Af in the formula, one is also free to choose S and 
Q, to depend upon e and x^(a, b, c, d) would still solve the difference equation. 

For example, let us consider 



X% = e'l^ 



AS 



X\ = e'l' 



X\ = 
so that 

{X\ - X%){X\ - K) - (A| - AD(A§ - Al) = e- 

(It is necessarily the case that (A4 — Af )(A3 — Al) is the difference of the other two.) 

Now, just for the sake of example, let us consider the solution which comes from 
using k — 2, n — A, 



/I \ 

1 e^^^ 

e^/" 

\0 e^/^ ei2A/ 



and 



foo\ 
1 
1 1 

\01J 
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This solution has the form^ 

x'^(a, b, c, d) — x\{a, b, c, d) — 2:2(0, 6, c, d) + 2:3(0, 6, c, d) 

where 

(a, 6, c, d) = e^/^ (l - e^/^) ' (l - e-''' - e^/^) ' (l + e^/- _ e^'') ' (l - e-^/^ + e^/- _ g^/-) ° 

xl{a, b, c, d)= (1- e^'^'Y (1 - e^/^ - e^/^)' (l + e'^' - e*^/^)' (l - e^/^ + gVe _ gS/e^ 

(1 - e^^/^)'' (1 - e^"/^ - ei5/^)' (l + e'^/' - e^'/'Y " + e'"/" - e^^/^)" 

and 

a:§(a, 6, c, d) = e"/^ (l - e«/^) (l - e^/^ - e«/^) ' (l + e^' - e«/^) ' (l - e^/^ + e^/^ - e«/^) " 

(1 - e^^/^)' (1 - ei°/^ - gisA)' (1 + ei4/^ - ei^/^)' (l - e'°/' + e^^/^ - e^s/^^ , 

Before even worrying about the structure of the equation, one must worry that 
this solution itself will not have a well defined or simple ultra-discrete limit. In 
particular, if we define 

u{a,b,c,d) = lim elnx^(a, 6, c, d), 

and we do not worry about the requirement that 13 7^ in Theorem 12. II then we 

would conclude that 

(11) 

u{a, 6, c, d) = Max(4+12a+126+12c-hl2d, 23a+235+23c+23d, ll+23a+236+23c+23d). 

(This is found by adding up for each factor in the expression of a;f with i = 1, 2, 3 
the highest numerator of the exponent when viewed as a fraction having e as the 
denominator. For instance, the first expression is found from x\ as 4 -I- 4c? -f 4c + 4 -|- 
4a -I- 8d -I- 8c -I- 86 -|- 8a.) Yet, as the examples of the previous section demonstrate, 
we ought to be worried that this formula is not accurate in the case that D = 
when Theorem 12. II is applied to this example. 

In fact, it turns out that there is no need to worry about this since it is always 
the case here that p| = 1. Note first that the third is always 11 more than the 
second, and so there is no possibility that these would be equal. By the same 
argument, there would be no problem if the first two were equal since they could 
not be the maximum (with the third being 11 more all of the time). Finally, we 
need only worry about the possibility that the first and third are equal. Certainly, 
it is possible for these two expressions to have the same value, but never for integer 
values of a, 6, c, and d. In particular, they are only equal when 

a = -|fe + c + d+^ 



^Up to a common multiple independent of a, b, c and d which docs not matter due to the 
bilinear nature of the equation. 
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A similar argument applies when we wish to compute the limit of the Hirota 
Bilinear Difference Equation to find the equation that this u(a, b, c, d) satisfies. 
Rewriting the HBDE as 

_ iu^a+l,b+l,c,d) + u^a,b,c+l,d+l)-4)/e , {u^a+l,b,c,d+l) + u^a,b+l,c+l,d)+eln{e^^^^ -e^*^^))/s 

— o-gg("°("+l'''>'=+l''') + "°('''''+l''=''^+l)+4)/e 

where u^{a,b,c,d) = e In |u(a, 6, c, and ai — ±1 as needed when u is negative, 
leads to its ultra-discretization as 

Max {u{a +l,b+l,c,d) + u{a, 6, c + 1, d + 1) - 4, u{a +l,b,c,d+l)+ u{a, 6 + 1, c + 1, rf) + 4) 

(12) = u{a + l,b,c+ l,d) + u{a, 5 + 1, c, d + 1) + 4. 

Then to show that (|lf |l gives a solution we apply Theorem 13.21 and note that since 
u(a + 1, 6 + 1, c, c?) = u(a, b, c + 1, d + 1) = u{a + l,b, c, d + I) = u{a, 6 + 1, c + 1, c?) 
the first argument of the Max operator is always 8 less than the second argument. 

5. Application: First Integrals 

In this section, we consider difference equations that admit first integrals and 
ask the following question: When does the ultra-discretization of the equation 
leave the ultra-discretization of the first integral invariant? The study performed 
in this section will be more algebraic than analytic: the results about first integrals 
will apply to the complete set of solutions as opposed to a particular solution as 
considered before. 

5.1. Preliminary examples. A well-known class of second order difference equa- 
tions is given by the QRT mappings discovered by Quispel, Roberts, and Thomp- 
son [14]. The ultra-discrete limit of members of that class of mappings together 
with their first integrals was first performed in [19] but, each time this is done, one 
has to check explicitly that the resulting ultra-discretization of the first integral is 
indeed invariant. 

Let us start with the following simple example of a difference equation 

/-,o\ aa;„ + l 

(13) Xn+l = , 

Xn^n—l 

where is constant in n. Equation is a member of the class of QRT mappings. 
As such, H13f) admits the following first integral 

a" 1 

(14) L {Xn-l,Xn) = Xn + Xn-1 H \ 1 . 

Xn Xn — 1 Xji Xn — l 

We say that t*^ is a first integral for (|f 3|l because i"^ (xn-i, Xn) = t^(x„,a;„+i) if x„ 
is a solution to (|f To obtain the ultra-discrete version of these two expressions, 
one writes a'^ = e^^^ and a;„ = e""/^. Then we take the limit as e ^ 0+ of the 
equation using the identity Q (note that in this particular example, we do not 
need the more general relation given by Theorem 12. 1|) . We obtain 

(15) Un+l = Max(u„ -I- ^, 0) - M„ - Un-l 

and 

(16) / = Max(M„, U„_i, A-Un-l, A~Un, -Un-Un-l). 
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The quantity / defined in H16|) is conserved by the equation H15|) . Following [9] , this 
can be shown directly by substituting the expression for w„_i in terms of ii„ and 
Un+l coming from (|15|) into the definition of / as follows: 

/„ = Max(^Un, -Un+l - Un + (m„ + A - U„, 

A + Un + l + Un - [Un + A)+, Un+l " (Wn + A)+^ 

= Max(^U„, A - Un+l, -Un - Un+l, A - Un, 

Max(^ + Un+l + Un - {Un + A)+, Un+l " + A) + )^ 

= Max{Un, A - Un+l, -Un - Un+l, A - Un, , Un+l) 
= Iji+li 

where we have introduced the notation {k)+ = Max (A:, 0). 

Another more complex example is given by a difference equation that is related 
(through the process of deautonomisation) to a discrete version of the third Painleve 
equation [6,15,16] 

(afa| + a;„) (af/a| + x„) 

(17) Xn+l - 



Xn-i (a\a\xn + 1) {a\la%Xn + 1) 



Again, as a member of the QRT family of mapping, (|17|1 has a first integral given 

by 

= afa^(l + af) {— h — ) a\a%{\^ Oif)(xn + Xn-i) 

e'i e € I - I 1 \,_e_e/ -^n , Xn — l 



+ a\ (X\a\ I XnXn-l ^ I + 0L-20l\ 

\ XnXn—1 / \Xn—l X 

The ultra-discretization of these two expressions is performed the same way as 
before by setting = e'^'/'^ and x„ = e""/*^. In the limit as e — > 0+, one gets 

(19) Un+l - 2Un + {Ai +A2- Un) + 

+ {Ai -A2- Un)+ - (Ai + A3 + Un) + 
-{Ai - A3+ Un)+ - Un-1, 

with first integral 

(20) I ^ Ma.x(^Ai + A3 - Un^i, Ai + A3 - Un, 

A1+A3 + 2A2 - Un-l, A1+A3+ 2A2 - Un, 
A1+A2+ Un, A1+A2+ Un-1, 
A1+A2+ 2A3 + Un, A1+A2+ 2A3 + Un-1, 
2Ai + A2 + A3 + Un + Un-1, 2^1 + A2 + A3 

Un Un — 1, 

A2+ A3 + Un - Un-1, A2 + A3 + Un-1 - M„ 



Checking that / defined above is a first integral is quite a cumbersome task. Instead, 
we will rely on a general theorem that we prove below. 
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5.2. Main result about first integrals. Consider the difference equation 

(21) ^n+d ~ fi^n- ^n+lJ ^n+d-l) 

where / is a rational function on M'* of the form 

(«m-=o'(a.r 

/(ao, ai, ad-i) = — ^ — 

A first integral for H21|l is defined to be a rational function l'^ on M"^ of the form 

(22) L^ao,ai,...,ad-i) ^ — -^-^^ — 

Efii (/3fato'(«.)'-) 

where hj, Uj G N, which remains invariant when restricted to any solution to (I21|l . 
that is, t^(ai, a2, a^-i, /(oo, Qi, Qd-i)) = '-'^(ao, ai, ad-i). Here, we assume 
the coefficients fi^ and /3| to be nonzero for e in a certain interval (0, (3). In all our 
examples, these coefficients will be exponentials and thus never zero. 

The following question arises: is the ultra-discretization of 1)22(1 a first integral 
for the ultra-discretization of (|21|l ? Precisely, does the function 

M ^1 A/ . 

/(6o,6i,...,5d_i) =Max(B, -I- > kjhj) -M&^{Bi + } kjbj), 

i—l ^ — ^ 2—1 ^ — ' 

where Bi , _B.; are defined as Bi = lim e In I /3f I and Bi ~ lim e In /3f , define a first 
integral for the equation 

(23) w„+d = g(u„, u„+i, Un+d-i) 
where 5 a function on R'^ defined by 

g(6o, 61, fed-i) = Max(Ai -I- > mijfej) - Max(Ai + > m^hj), 

i—l ^ — ^ i—l ^ — ' 

3=0 3=0 

where Ai are defined in Definition 1)3. 

Before stating the theorem, let us define the function D on R*^ by 

D{bo, bi, bd-i) = ^ CTic''' 

where 

P! (/(e''«/^e''l/^...,e''■'-l/^))'''-' 



(cTi is the sign of /3f (/(e''°/^ e''l/^ e^'^-i/'^))''''"' for e close to zero) and 
ajt = Ma.x{B, + } Ujbj+i + hd-ig{bo,bi, ...,fod-i)) 

i=l ' 

and 

d-2 

3 = + ^^ii^i+i + kd-ig{ba,bi, ...,bd-i) = M}. 
3=0 
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and 

eln\f3f\+j:''-^^k,bj+,+kd-ih%bo,h,...,bd-i)^Tl 
Ui = lim tor I £ J, 

e^0+ e 

where 



Note that 



lim If [bo, bi, bd-i) = g{bo, bi, bd-i). 



Furthermore, the function D is defined the same way by replacing a^, Hi^af, B^^ lij, M 
above by ct;, fit,af,Bt, kj , M 

Theorem 5.1. If D and D are nonzero almost everywhere in M."^, then I is a first 
integral for H^yp . 

Proof. By definition of /, we have the equality 



lim e In 



Furthermore, as a direct consequence of Theorem 12.11 since both D and D are 
nonzero almost everywhere in R'', the following equality must hold everywhere 
except in a subset of measure zero 

^lim £ln {^f (e^l/^ e^^/^ e'^-'/', ^h'iboM,-M-i)^^ ^ j^f^^^f,^^ ...^ ff(6o, 

Because both the argument of the limit and the right hand side of the equation are 
continuous in R'^, the equality holds in all R''. 

Since f is a first integral for (|21|l . the arguments of the two limits above are 
equal. Hence, since both limits exist, the right sides are equal which shows that / 
is a first integral for H23() . □ 

5.3. Examples. Consider the difference equation 

a" 



(24) Xn+2 = - Xn+l + 

Xn+1 

which admits the first integral given by 

(25) f = a;^+ix„ + Xn+ix'i - a^(x„+i + x„). 

(note that (|24|l was intentionally written so that Xn+2 is a function of x„ and Xn+i 
in order to apply Theorem 15 . II directlv with d = 2). Equation 124|l is related to a 
discrete version of the first Painleve equation [6,15,16]. For the ultra-discrete limit 
we choose = e^/^ and find 

(26) Un+2 = Max (u„, u„+i, A - u„+i) = g{un, u„+i). 



(27) /(u„, u„+i) = Max (2u„+i -I- m„, u„+i -|- 2u„, A + u„+i, A + u„) . 

The quantity defined in H27|l is not a first integral for 1)26(1 . To see this we con- 
sider the specific example with A < 8, in which we have that 7(3, 5) = 13 and 
7(5, g(5, 3)) = 7(5, 5) = 15 7^ 7(3, 5). 
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We now show that the conditions of Theorem 15. II are not met. Specifically, it is 
not true that the function D is nonzero almost everywhere in R^. We first set 

/5i = 1, ^1,0 — 2, ^1,1 — 1, 

P2 — li ^2,0 = li ^2,1 = 2, 
131 = ^3^0 = 1, ^3,1 = 0, 

PI = -e^/^ Z4.0 = 0, h,i = 1. 

and 

OT = Max (25(61, 62) + 62, 5(^1, 62) + 262, A + g{bi,b2), A + 62) • 
Furthermore, to apply Theorem 15.11 it is useful to write the following quantity 
explicitly 

Consider the case where 62 is greater than both 61 and A — bi. We obtain g(6i, 62) = 
62 and 3 = {1,2}. The sign of /(e^^^^^, e^^^^) for small positive e is negative which 
means that cti = — (T2 = 1. Since iii — 1^2 — 0, we finally get that D — e^^ —e^^ — 0. 
Thus the function D is zero in the semi-infinite triangular region of the plane 62 > &i 
and 62 > A — hi. The conditions of the Theorem l5.1l are thus not met. 

Let us now consider the two examples of Section l5Tl In both cases, all the cr^'s 
are equal to one. Thus, this automatically implies that both D and D are nonzero. 
In particular, this shows that (|20|l is indeed a first integral for (|19|l . 

We now consider a last example 

(28) Xn+2 -r—^^ Xn= f{Xn,Xn+l) 

with first integral 

(29) L = -l- a;„ — a;„_|_ia;„ — Xn+iXn- 
The ultra-discretization reads 

(30) u„+2 = Max (m„ -I- 2m„+i, m„+i, u„) - Max (2u„+i, 0) = 5i(u„, u„+i) 



(31) / = Max(2u„4.i, 2u„, 2u„+i + 2u„, u„+i + u„) . 

To apply Theorem l5.1l we set 

Z?! = I7 ^1,0 = 0, ^1,1 = 2, 

/3| = 1, ^2,0 = 2, ^2,1 = 0, 

/^a — —1, ^3,0 = ^34 = 2, 
/3| = —1, ^4^0 — ^44 = 1- 

and 

OT = Max (23(61,62), 262, 25(61, 62) + 262, 5(61, 62) + 62). 
Furthermore, it is useful to write the following explicitly 

> ~ I- p2b2/e 

To prove that the function D is nonzero almost everywhere in R^, we consider 
four cases. First, in the case 61,62 > 0, we have that 5(61,62) = 61 and U = {3} 
and thus, because there is only one term in Z?, it cannot be zero. In the case 
62 > 0, 61 < 0, 5(61,62) = Max(6i, —62) and J — {2}. Again, D cannot be zero. 
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In the case 62 < 0,6i > 62, we have that 17(61,62) = ^1 and 3 = {!}. Finally, in 
the case 61 < 62 < 0, 5(61,62) = 62 and 3 — {1,2,4}. Here, D is not trivially 
nonzero and a little more work is required. The sign of /(e''^^^, e^^/"^) for small e is 
positive, fii — 0, i = 1,2,4 and cti = 0-2 = ~o'4 — 1- Thus, = l + 1 — l = l7^0. 
We thus have proven that D is nonzero almost everywhere in which means, 
by Theorem 15 . II that H31|l is a first integral for (|3U|I . Note however that D is not 
nonzero everywhere in K^. For example, in the case 61 = and 62 > 0, we have 
that 3(61,62) ^0,3 = {2,3}, the sign of /(e''l/^ e''^/^) is negative, a2 = -0-3 = 1, 
and thus D = 0. 

6. Conclusions 

Although it is the case that problems can arise in computing a limit of the form 

lim e In T 

with T of the form ^ when some terms are negative, it is not true that these 
problems arise "most of the time" . Consequently, it is often still possible to use 
the procedure of ultra-discrete limits developed in [20] to produce integer solutions 
to equations involving the "Max" operator even in the absence of positivity, as 
illustrated in Sections El and 01 Furthermore, it is especially useful to apply the 
results of Theorem 12.11 to the question of which integrals of motion are preserved 
by an ultra-discrete limit since in that case one does not need to consider individual 
solutions (cf. Section 01. 

It should be noted that there are various different possibilities for the value of the 
limit in the case that D = and |3| > 1 which we do not presently have the ability 
to differentiate in general. For instance, in the example presented in Section 13.11 
the value of the limit is not the maximum of the exponents as usual but rather the 
next largest of the exponents. (This happens here because the sum that forms D 
does not only vanish in the limit but actually is equal to zero for all s, as briefly 
mentioned in [8].) On the other hand, in the case of the limit 

lim elnfe^i/^-e^^/^-He"^/^') 

with u\ = sin(2e)/£, U2 — 2sin(£)/e and u\ = \, the limit yields the maximum 
value of 2 despite the fact that Z? = 0. Most interestingly, in the example shown 
in Section 11.21 it is evident that the value of the limit is not determined by the 
values of the limits of the functions in the exponents, since the ultra-discrete limit, 
a — k, could be varied arbitrarily by selecting k without affecting the limits of the 
individual exponents, all equal to a. It would be interesting to extend Theorem l2.1l 
to be able to predict the value of the limit in the case that neither conditions (i) 
nor (ii) arc met. 

Acknowledgements: We are grateful to our colleagues Herb Silverman and Tom 
Kunkle for helpful discussions and advice. 
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